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A simple method for calculating the Casimir energy for a sphere is devel- 
oped which is based on a direct mode summation and counter integration in a 
complex plane of eigenfrequencies. The method uses only classical equations 
determining the eigenfrequencies of the quantum field under consideration. 
Efficiency of this approach is demonstrated by calculation of the Casimir en- 
ergy for a perfectly conducting spherical shell and fora massless scalar field 
obeying the Dirichlet and Neumann boundary conditions on sphere. The 
possibility of rationalizing the removal of divergences in this problem as a 
renormalization of both the energy and the radius of the sphere is discussed. 

12.20.-m, 12.20.Ds, 12.39.Ba 



I. INTRODUCTION 

The Casimir effect attracts the attention of investigators during last half of a century. 
More generally the Casimir effect can be defined as an influence of the boundness of the 
configuration space on the physical characteristics of the quantum field system (its energy, 
forces and momentum of forces acting on the boundaries and so on). 

This problem arises in different areas of theoretical physics: in quantum electrodynamics 
(attractive force between two uncharged conducting plates calculated by Casimir in 1948), in 
the theory of elementary particles (the bag models of hadrons treat the energy of quark and 
gluon fields located inside hadrons), in current cosmology, in physics of condensed matter 
(elucidation of the physical origin of sonoluminescence) and so on. 

When considering the Casimir effect different methods are used: Green function for- 
malism 1^, stress-tensor method 0, multiple scattering expansion 0, zeta regularization 
technique heat-kernel series 0, direct mode summation with counter integration 
Physical interpretations proposed in different approaches to calculating the Casimir effect 
are distinct. For example, the Casimir forces between two uncharged conducting plates can 
be treated both as the macroscopic manifestation of the van der Waals forces and as an effect 
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of zero point oscillations of vacuum electromagnetic field. In this situation it is worthwhile 
to separate in particular calculations an invariant, with respect to physical interpretation, 
"kernel" which gives the final result. 

In all the approaches to calculation of the Casimir effect a vague point is the procedure 
of unique separation and subsequent removal of the divergences. The lack of universal 
mathematically rigorous prescription for this purpose leads in some problems to different 
results when different methods are applied (for example, calculation of the Casimir energy for 
scalar massless field defined on the plane and subjected to the Dirichlet boundary conditions 
on a circle [Q). 

With regard to all this, the most simple, from mathematical point of view, methods of 
calculation of the Casimir effect have an obvious advantage because they right away allow 
one to reveal the difficulties generated by divergences. One of such methods is the direct 
summation of eigenfrequencies of quantum field system by making use of counter integration 
in complex frequency plane. For the first time this approach was proposed as a simple, in 
comparison with quantum field theory formalism, method of calculation of the van der Waals 
forces between dielectrics 0. Further it was widely used in other problems PJ^, p!0| , pn]1 . 



The main goal of this paper is to show the simplicity and efficiency of the direct mode 
summation by counter integration when calculating the Casimir energy for such a difficult 
boundary as sphere. This approach is completely based on using the classical frequencies 
of quantum field system concerned, and the main tool employed is the Cauchy theorem 
from complex analysis. In this approach the Casimir energy for perfectly conducting and 
infinitely thin spherical shell will be calculated. Then the Casimir energy of scalar massless 
field subjected to the Dirichlet or Neumann boundary conditions on sphere will be also 
derived. As far as we know the last problem (the Neumann boundary conditions) has not 
been considered in other approaches yet. Unlike other authors we propose to interpret the 
removal of divergences when calculating the Casimir effect for sphere as the renormalization 
not only of the energy but also of the radius of the sphere. 

The organization of the paper is as follows. In Sec. || we consider the vacuum energy 
of the electromagnetic field inside and outside the perfectly conducting spherical shell. In 
Sec. |IT1| the Casimir energy of massless scalar field obeying the Dirichlet or Neumann bound- 
ary conditions on sphere is calculated. Concluding remarks and discussion of necessity for 
renormalization of the sphere radius, in addition to the energy renormalization, are presented 
in Sec. IV, 



II. VACUUM ELECTROMAGNETIC FIELD INSIDE AND OUTSIDE THE 
PERFECTLY CONDUCTING SPHERICAL SHELL 

The starting point of our approach is the following definition of the Casimir energy 

E=Wi^s-^s). (2.1) 

^ s 

Here Us are the eigenfrequencies of the system under consideration, and lDs are those of 
the same system, when the parameters determining its boundaries take on some limiting 
values. We need equations for the oscillation frequencies of the electromagnetic field inside 
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and outside the perfectly conducting sphere with radius a. There are two modes of oscilla- 
tions: transverse-electric modes and transverse magnetic ones (TI5-modes and TM-modes, 
respectively). The eigenfrequencies of the TE-modes are defined by the equations [|12| 



ji{uja) = 0, 
h\^\uja) = 0, 

and the eigenfrequencies for the TM-modes are given by 

d 



dr 



d 
dr 



rhf'\ujr) 



0, 
0. 



(2.2) 
(2.3) 



(2.4) 
(2.5) 



In formulae ( |2.2| )- (|2!5| ) ji{z) and hf'\z) are the spherical Bessel functions [|T^ 



(2.6) 

and / = 1,2,.... Only positive roots of these equations Uni > 0, n = 1,2,... should 
be considered. Equations (|2.2|) and ( |2.4| ) specify the frequencies of the electromagnetic 
oscillations inside the sphere and Eqs. ( p.3|) and (|2.5|) give the frequencies outside the 
sphere 

In the case of spherical boundary the sum J2s ( P-l|) can be written as 



-| -| oo ( oo oo 

^E^^ = ^E E E^n^ = E(^ + l/2)5z, (2.7) 
^ s -^1=1 m=-ln=l 1=1 

Where Si = J2'^=i^nh and each frequency equation (^.2D-(P3D generates its partial sum 
Sr, a = l,...,4. 

For the partial sums S^"^ we use integral representation that follows from the Cauchy 
theorem flXll 



SI 



(a) 



— I dzz^\nf^"'\z,a). 
2Tci Jc dz 



Here f^°'\z,a) are the functions defining the frequency equations 



0, a = 1,2,3,4. 



(2.8) 

-( |2.5| ) in the form 

(2.9) 



The counter C encloses counterclockwise positive roots of these equations. Taking into 
account the position of the roots on real axis one can deform the counter C in such a way 
that it will consist of the segment [— ^A, iA] of the imaginary axis and a semicircle of radius A 
with A ^ cxD in the right half-plane. When A is fixed, the counter integral ( |2.8| ) gives the 
regularized value of corresponding frequency sum (it sums up the finite number of the roots 
of the frequency equation ( p.8| ) that lie inside the counter). 

From physical considerations it is clear that for negative values of the argument u the 
functions f^"\u, a) have to be defined by a condition 
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(2.10) 



This can be achieved, for example, by introducing a finite photon mass which is equated to 
zero at the end of the calculations. Separating the contributions of different parts of the 
counter C, we can rewrite formula (12. 81) as 



+A 



St^ = -^ I dy y ^ In f-\zy, a) 
2tt J ay 
-A " 



(2.11) 



Ca 



Here C\ is the semicircle of radius A introduced above. In the first term on the right-hand 
side of ( p. Ill) we can integrate by parts, the nonintegral term being omitted in view of ( p.lO[) . 
On the other hand this term can be removed by a subtraction which we shall discuss further. 

In accordance with the definition (|2.1|) in order to obtain a finite (observable) value of 
the Casimir energy it is necessary to perform the subtraction. As usual, we shall subtract 
the contribution of the Minkowski space that corresponds to the limit a = oo in Eq. ( p.llj) . 
Letting Si"'^ represent the value of the partial sum Sl"'^ which is to be subtracted from (|2.11|) 
we get 

+A 



St^ = ^ Jdy ln/(")(^2/,a^oo) 

-A 

+ — / 2rfln/(")(z,a-^ oo). 
27ri J 



(2.12) 



In view of an oscillating character of the function f^°'\z,a) (see below) we have on the 
semicircle Ca 



lini/(")(^,a) = /(")(^,a) 



(2.13) 



Thus, for the difference S^"^ — S}"'' in the definition of the Casimir energy ( p.lj ), we obtain 



=,(a) 



TX 



dy In 



f('^'){iy,a oo) 



(2.14) 



Here again the property (|2.10|) has been used. Now we proceed to substituting into Eq. (|2.14|) 
the concrete expressions for the functions /'•"^ defined by frequency equations ( p.2D -( ^^ . 
From Eq. ( p.2| ) we have 



/«(^2/, 



Ju{iya) 



Iu{ay) 



f^'^'){iy,a oo) lim Jy{iya) lim Iy{ay) 



27raye-''yiy{ay), 



(2.15) 
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where u = I + 1/2, and Iu{z) is the modified Bessel function Juiiz) = i^Ii,{z). We have used 
here the asymptotics of the function Iu{z) for fixed value of u and large z 

Uz) ^ (2.16) 



(2.17) 



From frequency equation ( p.3| ) it follows that 



a— >oo 



where Hj)^'> = Ju{z) + iNy{z) is the Hankel function of the first kind. Using the modified 
Bessel functions K^{z) = {7c/2y+^Hj^^'^ (iz) we rewrite Eq. (|2.17|) as 



my,a) Kiay) P-V.^a.^^^^^y (2.I8) 



/(^)(«?/, a ^ 00) \\m. K^jiay) V 



Here we employed the asymptotics of the function K^{z) for large z and fixed z/ |T3| 



KAz) ^ y^e-. (2.19) 



Thus, the TE-modes give the following contribution to Eq. ( p.l| ) 

^^^W _ ^{«)^ = -Jdy\n [2ayh{ay)K{ay)] . (2.20) 



«=i 



In the same way we deduce from the frequency equation 

f^^\iy,a) Jy{iya)/2 + iyaJl{iya) 



f^^\iy,a^ 00) limJ^J^iya) /2 + iyaJ^iya)] 

^ Iuiya)/2 + yaIl{ya) 
lim [lAya)/2 + yarAya)] 



(2.21) 



The prime over the Bessel function means the differentiation with respect to its argument. 
From (|2.1tj|) it follows that 



lim [lAay)/2 + aylliay)] = Me^y. (2.22) 



In view of this, formula ( p.21| ) assumes the form 



/(3)(^y,a) _ IEe-^y[I,(ay)/2 + ayIliay)]. (2.23) 



/(3)(i?/,a ^ 00) \ ay 
In the same way, for the frequency equation ( |2.5| ) we obtain 
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/{4)(2y,a^oo) 

Hi^\zya)/2 + iayH(^y{zay) 



lim 



Hi^\iay)/2 + iayHl'' (iay) 



Ku{ay)/2 + ayKl{ay) 
lim [K^{ay)/2 + ayKl{ay)] 



Tiya 



-e'^y [K,{ay)/2 + ayKl{ay)]. 



(2.24) 



Summing up Eq. ( p.23| ) and ( |2.24| ), we arrive at the contribution of the TM- modes 

4 



«=3 



(a) n(a); 



1 
TT 

X 



j dy\n\- 



2 

ay 



-Iy{ay) + aylKay) 



-Ku{ay) + ayKKay) 



Finally for the Casimir energy (|2.1|) we obtain from (|2.7| ), (|2.2(]| ) and (|2.25 ) 



(2.25) 



E 



na 



1=1 



Y.[^ + ^)Jdyln [l-{a[{y)y 



where the notation 



-I,{y) + yll{y) 



(2.26) 



is introduced. Using the value of the Wronskian of the modified Bessel functions Iu{y) and 
KM 

iMK'M - I'MKM = -- 

y 

one can show that 

ai{y) = yI,{y)Kuiy), p = 1 + 1/2. 

The integral in ( p.26|) converges. This follows from the asymptotics of (j[{y) for large y and 
fixed z/ = / + 1/2 



1 



2y^ 



1 4z/2 - 1 

,^"2~(2^ + --; 



(2.27) 



Formula ( |2.26| ) coincides with Eq. (5.1) in paper |1[|, on the condition that the cut-off factor 
in the last equation is omitted and the integration by parts is performed. To carry out the 
summation with respect to / in ( |2.2(j| ) one needs the behavior of the integral 



6 



Qi 



1 + 1/2 



TT 



dyln 1 - {(y[{y)f 



(2.28) 



at large /. Applying the uniform with respect to z asymptotics for the modified Bessel 
functions at large z/ [jl|,|l3| 



h{uz)K^{uz) ^ Y.JlTzW' 



(2.29) 



we obtain from ( p.28| ) 



Q CXJ 

u f 
Qi ^ — dz \n 

TT J 



4z/2(l-22) 



dz 



Att J (1 + z 



2\3 



, / —>■ OO. 

64' 



(2.30) 



Thus, the sum ( p.26|) at large / diverges as + 1/2)°- To determine the finite value for 

this sum we rewrite ( |2.26| ) in the following way 



-| OO 

E = -y: 



OO 

a 



■ 3 3 



where 

Qi = Qi + 3/64. 
The sum Qi converges because at large I 
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(2.31) 



(2.32) 



Qi 



16384Z/2 

The last divergent sum in (|2.31|) can be defined by using the Hurwitz zeta function |14 



2.33) 



OO -| 

c(^,g) = E 



(? + '^)^' 



which at g = 1/2 is related with the Riemann ^-function |]T5 

C(z,l/2) = (2^-l)C(z). 

From (|2]3|) it follows that 



z=i 



64a 



^ (C(0,l/2)-l) ^ 



64a 



(2.34) 



(2.35) 



(2.36) 
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since C(0, 1/2) = 0. 

Finally for the Casimir energy we obtain 



1 OO _ o 

E = -T.Qi + ^^ (2-37) 

where Qi is defined in ( p.32|) and (|2.28|) . The sum Qi iii the right-hand side of ( |2.37| ) can 
be estimated with allowance for the asymptotics (|2.33|) 

|:e.--i^i:(j^ = -^ic(2.i/2)-4] 

^-Jlfy -4) =-0.000514.... (2.38) 

Thus, the main contribution to ( |2.37| ) is given by the second term and to a good approxi- 
mation one can put for the Casimir energy 

E^ — = - 0.046875. (2.39) 

64a a ^ ' 

Taking into account ( p.38| ) we get instead of ( p.39|) 



E~ -0.046361.... 

a 



With greater accuracy this energy has been calculated in |]1|. 



III. SCALAR FIELD OBEYING THE DIRICHLET OR NEUMANN BOUNDARY 

CONDITIONS ON SPHERE 

The suggested method can be easily applied to the calculation of the Casimir energy of 
a massless scalar field subjected to boundary conditions on a sphere. Let us first consider 
the Dirichlet boundary conditions. In this case the eigenfrequencies are given by Eqs. ( p.2|) 
and dU) with / = 0, 1, 2 . . .. From (|2l|), (|3) and ([OOl ) we obtain 



E(^) = Y^[i + -)^-j dy \n[2ayh{ay)K,{ay)] 



1=0 

00 





00 



1 / 1\ 1 /■ 
= -T.{^+o)- dy \n[2yI,{y)KM], (3-1) 

u = l+ 1/2. 

It is worth comparing this expression with Eq. (3.5), derived in paper [0] by the Green 
function method. The essential advantage of our approaches is absence in the integrand of 
the constant terms which lead to the divergences (so-called contact terms). For large u the 
integral in (|3.1|) behaves as follows 
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Qi ^ I dy \n[2yhiy)KM] 



,,2 



+ 



35 



128 32768Z/2 



(3.2) 



The first two terms in ( |3.2| ) give rise to divergences when summing up with respect to I in 
These divergences can be removed as in the previous section by means of the Hurwitz 
C-function (|2.34]) . On this purpose we rewrite (|3.1|) 



-| oo 
^ 1=0 

1 /, 1 
> / + - 



128 



128a V 2. 



(3.3) 



Here we have added and subtracted under the sum sign two first terms of asymptotic ex- 
pansion ( p. 21) . Taking into account (|2.34|) we obtain 



EP) = lVg,_lc (-2,- 
af^^^ 2a \ '2 

-| oo 
" 1=0 



128a V 2 



(3.4) 



where 



Qi 



+ 



128 



(3.5) 



0. 



We derived the last equality in ( |3.4| ) bearing in mind that C(~2, 1/2) = and C(0, 1/2) 
By virtue of the asymptotics (|3^ ) the last sum in (|3.4|) obviously converges. 

With increase of z/, the modified Bessel functions Iu{,z) and Ky{z) are rapidly approaching 
their uniform asymptotics. That is why even for comparatively small values of / we can 
assume with allowance for (|3.2|) and (|3.5| 



asym 



35 



32768i/2 



(3.6) 



This simplifies the numerical calculations considerably. It is important to note that direct 
calculation of the modified Bessel functions Iu{z) and K^{z) with desired accuracy for all z 



and at large u is a technically difficult problem ||T^. It is especially concerns the product 
I^{z)K^{z) encountered in the integrand of (|3.2|) . 

Table | demonstrates the applicability of the formula ( |3.6| ). Calculating Qi for / < 3 
by means of numerical integration of the expressions including the product Iu{y)Ki_,{y) (see 
Eqs. ( |3.5| ) and ( |3.2| )) and using the asymptotic formula ( |3.6| ) for / > 3 we obtain 
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E(^) = 1 0.002819. . .. 

a 



(3.7) 



With greater accuracy the Casimir energy for spherical conducting sliell was calculated 
16| by making use of the Green function technique. 



m 



Now we proceed to the consideration of the energy of a scalar field obeying the Neumann 
boundary conditions on the sphere. The eigenfrequencies inside and outside the sphere are 
defined now by equations 



^ [3l{ujr)] \r=a = 0, 



d_ 

dr 



hf\ur) 



(3.8) 
(3.9) 



where / takes values 0, 1, 2, 



By analogy with Eq. (|2.25|) we derive 



°° f l\ I [ ( 2 rl , 

5^ r + o ~ / dy\n\ i^Iuiay) - aylUay) 

f^o^ 27 TT 7 [ ay [2 



X 



-K{ay) - ayKl{ay) 



(3.10) 



As in the case of electromagnetic field it is convenient at first to consider the sum 
From (13) and (|310|) it follows that 



1 °° / ^ \ °r r 



1=0 



(3.11) 



where the notation 



1 



-JAy) - yKiy) 



KM-yK'M 



(3.12) 



is introduced. It is easy to show that 



d (I 



(3.13) 



The integral in ( |3.11| ) converges since for large y and fixed / the following estimation 



^^l{y) - y ^oo 
y 



(3.14) 



holds. 

The convergence of the sum in ( |3.11|) is determined by the behavior at large / of the 
expression 
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Pi = ^-^^ I dy\n . (3.15) 



TT 





Using the uniform asymptotics of the modified Bessel functions [13|, we obtain 



Pi^--- + 0{iy--'). (3.16) 

64 163841/2 ^ ^ ^ ' 



The first term in this expression leads to a divergence in ( p. Ill) when summing with 



respect to /. We again overcome this difficulty using the Hurwitz zeta function 

£(« + £(A^),if;p, = if;p,_i?_cfo,i 

/=( 

oo 



«/=o ««=o 64a V 2 



where Pi is the renormahzed value of Pi 

19 

P^ = P^ + ^^ (3-18) 

When calculating the last sum in (|3.17|) we can again assume 

1 53 

^'-^r™ = -Y^^> - = 1 + 1/2. (3.19) 

Numerical calculation shows (see Table ^ that with increasing I Pi approaches rapidly 
to its asymptotic value defined by Eq. ( 3.19|) . If the required accuracy is not too high one 



can use Eq. ( ^.19| ) even for / > 3. This leads to the result 

+ = _i 0.220958... . (3.20) 



Hence, taking into account ( |3.7| ) we have for the Casimir energy of massless scalar field 
obeying the Neumann boundary conditions on the sphere 

^(A^) = _1 0.223777. . . (3.21) 
a 

As far as we know this result is obtained here for the first time. 



IV. CONCLUSION 

The direct summation of eigenfrequencies by calculation of the Casimir effect for nonflat 
boundaries (specifically for sphere) has been used only in pioneer paper by Boyer [^]. The 



fact that done by us is actually a development and maximum simplification of the Boyer 
method and bringing it to such a form when numerical calculations are practically not re- 
quired (see Eq. (|2.39| )), and, what is more important, cut-off functions are not used. In other 
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approaches, for example, by making use of the Green function formahsm [ll|,[T6|, transition 
to the imaginary frequencies is used without detailed justification. Contour integration in 
the mode summation method supplies a clear explanation for this technical trick. 

Let us turn now to removing the divergences in the problem under consideration. Once 
renormalization of the Casimir energy by the formula ( p.l| ) is accomplished the divergence 
for all that remains. In the general case the formula concerned has the form 



E = ^ ^, 4.1 

a 

where Ci is finite constant, and C2 is a divergent expression. Thus for example, for the 
electromagnetic field Ci = 0.046176 . . . and is given by the divergent series 

q 00 

C2~ = -7^E(^ + l/2)°. (4.2) 

To remove this divergence we have applied the formal technique of the zeta function renor- 
malization. In Ref. 0] the finite result was obtained here by making use of an exponential 
cutting multiplier splitting the arguments of the field operators in energy-momentum tensor. 
When calculating the Casimir energy for scalar massless field obeying the Dirichlet boundary 
conditions on sphere Ci = 0.002819 . . . and stands for the sum of the divergent series 

1 00 -1 00 

= -o + 1/2)^ - + 1/2)°. (4.3) 

2 1=0 128 

Both in our paper and in Ref. this divergence has been taken away by formal technique 
of the Hurwitz (^-function. 

Here the following question arises: to renormalization of which parameter does the re- 
moval of the divergence correspond? After renormalization of the Casimir energy accord- 
ing to Eq. (|2.1|) only one parameter, namely, radius of the sphere a, is available. Therefore 
it is natural to treat the removal of the divergence as the renormalization of the sphere 
radius. This can be done in a standard way by transition from the initial (bare) radius 
a to the physical (observable) radius aphys- a = CLphys + ^cl, where Sa is the appropriate 
counterterm. In view of this, Eq. ( |4.1| ) can be rewritten as 

O O-phys + ^0- O-phys (1 + ^O'/O'phys) 



Setting 



5a C2 



(^phys C\ 



00 

5 



(4.5) 



one arrives at the finite result 



E=-^. (4.6) 

^phys 
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This reasoning seems to be more consistent as compared with, for example, introduction 
into consideration of a phenomenological interaction ||^ locahzed on sphere with subsequent 
renormahzation of the couphng constant of this interaction that finally takes up divergence 
C^. It stands no reason that explanation suggested is also applicable to the calculation of 
the Casimir effect for the field confined inside the cavity. Certainly, in this case the values 
of the constants Ci and in Eq. ( [4.1|) will be different. 
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TABLE I. Qi obtained by numerical integration according to Eqs. (3.2), (3.5), and asymptotic 
formula (3.6) for Q^^^"^. 



I Qi grz 

0.001913 0.004273 

1 0.000398 0.000474 

2 0.000159 0.000171 

3 0.000084 0.000087 

4 0.000052 0.000053 



TABLE II. 


With increasing I Pi rapidly approaches 


to its asymptotic value p^^^^ given in 


(3.19). 






I 


Pi 


pasym 





-0.211491 


-0.037351 


1 


-0.004800 


-0.004150 


2 


-0.001582 


-0.001494 


3 


-0.000775 


-0.000762 
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